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A class of generalized non-minimal coupling theories is investigated, in search of scaling 
attractors able to provide an accelerated expansion at the present time. Solutions are found 
in the strong coupling regime and when the coupling function and the potential verify a 
simple relation. In such cases, which include power law and exponential functions, the 
. . . dynamics is independent of the exact form of the coupling and the potential. The constraint 

0^ \ from the time variability of G, however, limits the fraction of energy in the scalar field to 

0^ . less than 4% of the total energy density, and excludes accelerated solutions at the present. 
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I. INTRODUCTION 



A great deal of effort has been devoted in recent times to the dynamics of scalar fields in the radiation 
^ ' and matter dominated era. The motivations are manyfold: first, several theories of fundamental physics 

, predict the existence of scalar fields [0 [§) second, a slowly rolling scalar field may mimic the behavior 

of a cosmological constant at the present time, in agreement with apopular model of structure formation 
\ and with the observation of an accelerated space expansion [g ; third, the scalar field may alleviate 

■ the constraints on a true cosmological constant Q ; fourth, the additional source of fluctuations produced by 

' the scalar field may give new observable effects on the cosmic microwave background and on the structure 

g^ : formation % % %. ______ 

So far, most work focused on fields with minimal coupling to gravity W (s) M H] |10[. In this case, 
the dynamics in a homogeneous and isotropic space-time is completely determined once one specifies the 
matter fluid equation of state and the field potential. For the former, the obvious choices of interest are 
O \ the equation of state of a relativistic fluid and of a pressureless one. For the latter, although there are no 

' observations or fundamental principles to guide our investigation, potentials like power-laws, exponential and 

^ \ a handful of other cases have been selected, basing either on simplicity or on some particle physics model. 

Among the infinite solutions of the system of equations, the attractor solutions are of course of the greatest 
interest. Among the attractors, those which have a power-law behavior, denoted also as scaling solutions, are 
particularly simple to find and to study. Consequently, the study of the scalar field dynamics has focussed 
\^ ' on the search of scaling attractors. To be interesting for cosmological purposes, these attractor solutions 

5^ , must also lead to a energy density in the scalar field which is a non-negligible fraction of the total energy 

density. Finally, if we want to explain recent observations of the large-scale geometry of the space-time |Q 
p^ , the scalar factor has to be accelerated at the present. 

In minimal coupling theories the Lagrangian is the sum of the Einstein-Hilbert gravity Lagrangian and of 
the scalar field sector. The non-minimal coupling (NMC) adds a new term which, in its simplest form, may 
be written as (for a more general form that includes derivatives see Ref. JTsf) 

fim- (1) 

For instance, Refs. H |1 @ adopted / ~ 0^ discussing the model in the context of infiation. In 
Ref. several attractor scaling solutions in the matter dominated regime with power-law and exponential 
potential have been found. Other forms of /((/)) have been considered (see e.g. ]l9|). 

A common feature of all these investigations, perhaps obviously, is the choice of specific potentials and 
coupling functions. The purpose of this paper is to show that it is possible to find attractor solutions in NMC 
models in which both the coupling f{(j)) and the potential V{(l)) are left unspecified, and only their relation 
matters. In other words, we will find a class of models in which the dynamics of the system is independent 
of the coupling and of the potential, and depends only on their relation. In particular, we will find attractor 
solutions for all models for which wc can write |Q 

V{4>)^Af{4>Y'. (2) 

This relation holds, for instance, when both V and / are power-law, or exponential, but is also valid for much 
more complicated functions, like products of power law and exponential. In the limit of strong coupling, 
the dynamics of the cosmological solution will be shown to depend essentially only on M and on the fiuid 
matter equation of state. 



After performing a conformal rcscaling of the metric, the NMC system is written as a scalar field in pm'e 
General Relativity with an exponential potential and an extra coupling to the ordinary matter (see e.g. 
pH l , in which however only the case / ~ (/>^ has been discussed). This system shows a surprisingly rich 
phase space structure, with four different attractors. Two of these are qualitatively similar to the attractors 
found in the system without extra coupling. The other two however are new, and have not been previously 
identified. Although we derived this system from a class of NMC theories, we remark that it is interesting 
on its own, and many cosmological properties of its trajectories have yet to be worked out. Here we study 
it mainly to constrain the NMC model, and find that the constraint on the variability of the gravitational 
constant rules out this class of models as explanation for the accelerated expansion rate of the Universe. 

In the next section we work out the field equations. In Section 3 we find and discuss the attractor solutions, 
in Section 4 we discuss their cosmological properties, and in the final Section we draw the conclusion and 
point to new developments. 

II. FIELD EQUATIONS 

Consider the Lagrangian of a NMC scalar field plus a perfect fluid matter component (k^ = SttA/j^^) 

Ltot = L{(t}, R) + 2k^L^ + 2K'^Lmatter, (3) 

L{^,R) = ~fi^)R, (4) 
= l^,^.^'^ ~ Vi^). (5) 

Contrary to the usual notation, we found convenient to include into the constant that produces the 
Einstein-Hilbert term, so that our / is 1 + in the notation of, e.g., Ref. ||l6| jl^l). We will 

always assume / > 0, since it acts as an effective gravitational constant, 

Ge// = {2^^f)-'. (6) 

The Einstein equations are 
■f 



7:9iJ.v{L — L rR) — g^u^L^n + {L^b);^v + ^ + 



(7) 



where L r denotes here dL/dR, and where the scalar field energy-momentum tensor is 



and the fluid tensor is 

Tf,u{ni) = {p+ p)u^Uy - (9) 
Now, under the conformal transformation 

g^iv = e^"5p^, (10) 
the following transformations (see e .g. [^) occur: the kinetic term 

K^„y(<p) = (t),^,(t).M - ^5'm'j'/',q0'", (11) 

remains invaried (-ftT^^ = K^^); the potential term gf^vVi^') becomes e~'^'^g^yV{(l)); and the perfect fluid 
tensor becomes 

^^uirn) ^ ^^firn)- (1^) 

Putting 

2a; = log/, (13) 

it follows that the equations in the rescaled metric (sometimes called Einstein frame, while the old metric is 
the Jordan frame) are 



(14) 



where 



F\<l>) 



(15) 



where = 2k^/3 and where the prime denotes derivation with respect to 0. We can then define a new 
canonical field 



I dm), 



a new potential 



and a new matter tensor 



Finally, all these definitions lead to the canonical equations in the new metric g^n 

The new matter energy-momentum tensor can be written as 

f;(;„) = diagipe-'- -pe-'"^ , -pe~'- , -pe'^^ = diag{p\ -p* , -p* , -p*). 
As a last step, we rewrite the new metric in the Fricdmannian form 

gf^„ = diag{l, -a^ - a^, -a^), 
where the old time and the old scale factor are 

t = j e-'^^^Ut, 

and 



(16) 

(17) 

(18) 

(19) 

(20) 
(21) 

(22) 
(23) 



The equation of motion for the fields are obtained as the covariant conservation laws of the energy tensors. 
In the old frame they read 



The transformation to the new frame is performed according to the rules 

R = e^'^iR - 6g"^uj^aLo,(i + 6Da;), 
□ = e2'^(a-2g"VaV3). 



(24) 



(25) 



From now on, we omit all the tilde, until we return to the original quantities. Finally, the full set of equations 
in the Priedmann metric read: 



p*+3H{p*+p*) =-lw,^4,{p*-3p*), 

3H^ = K^(p* + iy? + u) , 



(26) 
(27) 

(28) 
(29) 



where 



iy = log /((/.), W,^, = ^. (30) 

As already remarked in the Introduction, the system ( p6[ - |29| ), here derived from a NMC model, is interesting 
on its own. Indeed, we can regard either the Jordan or the Einstein frame as the physical one. In the 
former case, we have to express the solutions of the above system back in the original frame, and study 
its cosmological consequences in the original frame, as we will do below. In the latter case, the solutions 
of the system are the physical solutions, and their properties can be directly compared to observations. In 
particular, the constraints from the variability of G , which we will find to limit heavily the cosmological 
viability of our solutions, apply only assuming the physical frame to be the original Jordan one. 



III. SOLUTIONS 

The full dynamics of the system ([2^p9|) is specified by the potential U and by the equation of state 
p = {w — l)p. In the following we consider only < w < 2. To write down the potential U{'tp), we have 
first to find the relation between ip and (j). This is where the possibility of a dynamics independent of the 
potential and of the coupling function arises. In fact, if we assume that pO| 

/" » cV, (31) 

then we can simplify Eq. ([l5|): 

^'(0)=(^^j • ^^^^ 

It follows 

c^ = J J=^''&f^ (33) 

where the integration constant can be absorbed into a redefinition of ?/; . It follows that the conformal 
function uj equals cip/l. Therefore, once we have the dynamics of if) in the transformed metric, we can write 
down the solution in terms of the original metric without having to specify /(</>), provided we express also 
the potential V{cj)) as a function of f{(j>). With the assumption ( |3l] ) we get W^^ = c in the system ( p6[ - p9| ), 
so that putting 

/? = 4 - 3w, 

the first two equations become 

p* + 3Hwp* = -^c(3tPp*. (34) 

The condition (|3^) holds true in several cases. For instance, it is verified for large (f> by any function /(</>) 
which grows faster than quadratically, that is lim^^ao f {(j)) / 4''^ ^ oo. In the often-studied quadratic case, 
/ = 1 + for large (j) we can put f((j)) = . Then, instead of Eq. (|3ll) , one has /'^ = 4k^^/ 



and all that changes is that in Eq. (32) and Eq. ( [33| ) is replaced by c^/(l + 1/6^). In this case, all 
the results found below become exact. The weak coupling limit in which K'^S^cfP ^ 1, i.e. <C (k^c/P)^^, 
on the other hand, is excluded in the present analysis. We could then label our case as the strong coupling 
limit. In fact, it is easily seen that it corresponds to the limit in which the Lagrangian can be approximated 
as —f{(t))R — 2K^V{(j)), neglecting the kinetic term ■ Notice however that this does not imply that 

the scalar field kinetic terms in the field equations are negligible, because the non-minimal coupling itself 
introduces other kinetic terms. 

Now, as anticipated, suppose we can write V{4)) — Af{(j))^^ .The potential becomes then 

UW^^^M^=AeVv'^^^^, (35) 



where 



fi = M -2. 



(36) 



Therefore, the potential can be written as an exponential, whatever the shape of V and of /, provided 
that the condition ( pT| ) and the relation (||) are fulfilled. The sign of fj, selects the direction in which the 
field ■0, and thus the variable /, rolls. U fi > 0, tp rolls toward — oo, so that / ^ 0, and the effective 
gravitational constant Ge// increases with time. In the opposite case, /it < 0, we have that Geff decreases in 
the future. We emphasize that if /(0) is quadratic, then all results below remain valid provided /i is replaced 
by iig = + 1/60^/2 and /? by = /?/(l + \/&(,Y'\ 

The scalar field dynamics in NMC theories is then reduced to the scalar field dynamics in pure general 
relativity with an exponential potential and with a scalar field/matter coupling. In the radiation case in 
which w = 4/3, we have /? = 0, and the source terms decouple. The decoupling occurs also when we can 
neglect the matter energy density Vlp- = K?p* /?>H'^ with respect to the scalar field energy density 11^. In 
these cases, the problem is identical to that already solved in, e.g., Ref. Q |9| The case /3 7^ 

has been already discussed by Wetterich in where some of its attractors have been identified. Here we 
extend the analysis to the full classification of critical points and attractors (finding two new attractors) 
and express the solutions in terms of the old frame. We keep (3 as an independent parameter as long as 
possible, and proceed to replacing it by 4 — iw only in the graphics, in order to narrow the parameter space 
to two dimensions, namely w and /i. The formulas apply however to the more general case, unless otherwise 
specified. 

Following Copeland et al. IM we define 



Kip 

7^' 



(37) 



and introduce the independent variable a — \oga(t). Notice that and give the fraction of total energy 
density carried by the scalar field kinetic and potential energy, respectively . Then, we can rewrite the 
system (26 -291) as 



x' = —3a; + 3x 
y' = fixy + 3y 



x^ + ^w{l ~ x^ - y^) 
x^ + ^w{l ~ x^ - y^) 



t^y^ + ^m-x^-y^), 



(38) 



where the prime is here d/ da. The system is invariant under the change of sign of y and of a. Since it is 



also limited by the condition p* > to the circle x 
positive y. The critical points, those that verify x' ~ 
equation of state is 



+ y < 1, we may study only the unitary semicircle of 
y' = 0, are scaling solutions, on which the scalar field 



2a;2 

= ^ = const, 

x^ + y^ 



(39) 



the scalar field total energy density is fi^ = x + y , and the scale factor is 



2 

3w 



w 



w + U,^{wtp — w) 



(40) 



(the slope p is not to be confused with the pressure) . 

Copeland et al. |0 have shown that the system (^8|) with /3 = and an exponential potential has up to 
five critical points, that can be classified according to the dominant energy density: one dominated by the 
scalar field total energy density (let us label this point as solution a and refer to its coordinates as a;a,ya ), 
one in which the fractions of energy density in the matter and in the field are both non-zero (labelled 6), 
one dominated by the matter field (c), and finally two dominated by the kinetic energy of the scalar field, of 
which one at a; = — 1 {d) and one at a; = +1 (e). 

The critical points on which the matter field becomes negligible reduce to the /? = case: therefore, the 
solutions a, d, and e remain the same also for /9 ^ 0. The points h and c are instead modified. The solution 
c is no longer matter dominated: rather, the scalar field kinetic energy and the matter energy take up a 
constant fraction of the total energy. In MDE, the scalar field kinetic energy amounts to fi^ = 1/9. The 
critical points in the general case /? 7^ are listed in Tab. I, where we put g{l3, w, /i) = + 2/3/i + 18w. 
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Tab. I 



Although the number and position of the critical points is affected only quantitatively by the extra coupling, 

their stability properties arc modified in a more radical way. In particular, while for /? = only the points 
a and b can be attractors, here we show that also c and d may be stable. Only the point e remains always 
unstable. 

The stability analysis is performed as usual by linearization around the critical points. The parametric 
regions in which the real part of both eigenvalues of the linearization matrix is negative are regions of 
stability. To simplify the discussion, we only consider the crucial property of stability versus instability, 
paying no attention to the topography of the critical point (whether it is a knot, spiral, or saddle). In the 
following, we say that an attractor exists if it lies in the region Q < + < 1. The parameter spaces are 
plotted in Fig. 1. 

Point a . 

The solution a exists for < 9, and is an attractor only for < fj, < fj,+ where 



(e.g., /x_ = —2.39 and /i+ = 1.89 for w = 1). On this attractor we have = 2/^^/9 and 

3/m', 



Pa 



(41) 



(42) 



inflationary if |/Lt| < VS- 
Point b . 

The attractor b exists and is stable in the region delimited hy fi < fj,- and IJ, > IJ.+ and the two branches 
of the curve 



1 



Mo = -^(/3' + 18«;-9«;2). 



The scale factor slope on the attractor is 



Pb 



^(1 + A 
3u; V 2/i 



and, for /? = 4 — 3w, is inflationary within the two branches of the curve 

4-3w 

Mi 



3u;- 2 



(43) 



(44) 



(45) 



It is remarkable that the inflationary region for the point b includes values smaller than w w 0.91, and 
therefore excludes the MDE equation of state w = 1. This conclusion is not changed by replacing /x and /? 
with their counterparts fXg and 0q in the case of a quadratic coupling f{(p). 
Point c . 

This point exists for w < 5/3, and is stable below the lower branch and above the upper branch of /kq . 
The slope is 



Pc 



6(2 - w) 



/32 + 9(2-«;)w' 



and it is never accelerated if /3 = 4 — 3w. The point c shares with b the property that matter and scalar field 
have both a non- vanishing fraction of the energy density. 
Point d . 

This point exists for all values of the parameters, and if /? = 4 — 3w is stable for w > 5/3 and > 3. Its 
slope is always pa = 1/3. 
Point e. 

This point exists and is unstable for all values of the parameters if /? = 4 — 3w. 



The complex structure of the parameter space is summarized in Fig. 2. Notice that 1), for each value 
of the parameters w,^ there is one and only one attractor; 2) for w = 1 the points a, 6 or c can be stable, 
depending on fi; 3) these solutions are inflationary in the shaded region; 4) only the point a can be accelerated 
for w; = 1 or larger. In Fig. 3 we present four phase spaces displaying in turn the four possible attractors. 
The parameters correspond to the points marked with stars in Fig. 2. As already remarked, attractors c and 
d have not been previously noticed. Also, it is important to remark that the attractors are not only locally 
stable, but extend their basin of attraction to all of the phase space. That is, any possible initial condition 
lead to the attractor. 



IV. BACK TO THE JORDAN FRAME 



Here we leave the dynamical analysis of the system in the rescaled frame and get back to the original one. 
What the attractors look like in the Jordan frame? 

Reintroducing the tildes, we have along the attractors a and b (from now on, quantities without tildes are 
expressed in the original metric) 



(for /i ^ 0) where 



ci; = --\og\t/Ta,b\, (46) 



r-l = V2A^,c^. (47) 

Va.b 



On the attractors c and d , for which y — 0, ip —^ oo, and the conformal transformation cannot be performed. 
Since ijj is proportional to log/, the attractors c and d lead to an effective gravitational constant that is 
either zero or infinite and are therefore to be rejected as possible solutions in the Jordan frame. Of course, 
trajectories that have not already reached the attractor cannot be excluded, but these are not scaling 
solutions, and will not be further considered in this paper. 
Form Eq. ( ^6| ) it follows (neglecting the subscripts) 

e2'^ = (f/T)-^^. (48) 

From the latter expression we can evaluate the relation between the old and new time and scale factor, given 
by Eq. (|2|) and (|| ). We obtain (for /i ^ 0, -1) 

— fi 

a^i-i/(i+M)a(i). (49) 

As can be seen, for /i ±oo the old and new metric coincide; in this limit the scalar field vanishes on the 
attractor, and the system reduces to the pure perfect fluid Friedmann case. 
It follows that in the original variables the scale factor is again a power law 

a ~ , P = — . 50) 

1 + M 

On the attractor a, pa — i/ p^ , which is inflationary (both in the original and in the rescaled frame) if /i^ < 3, 
that is 2 - V3 < M < 2 + V3 , we have 

/ /^ + 3 , . 

Pa = /-I , V 51 

m(1 + N 

Consider now some special cases. If > /i > — 1 the scale factor follows a pole-like inflation, a ^ (fp — i)^° 
with negative exponent. For = ,(i.e. V ~ p) the old and new metric coincide (up to a constant), the 
field freezes to a constant and its energy drives a deSitter expansion. The system reduces asymptotically to 
pure general relativity with a cosmological constant. Finally, for /i = — 1 (i.e. V ^ /), the scale factor is 
power law accelerated in the new frame, but maps again to a deSitter expansion in the original frame. If / 
is quadratic in then the inflationary condition on the solution a reads 

< 3(1 + 1/60- (52) 



On the attractor b , on the other hand, putting pi, — 2/ {3w') with w' = w {1 + ^ we obtain 



pj, — » 2/3w for fi ±00 , as expected. Since the property of being accelerated is conformally 
invariant (for positive definite conformal factors), going back to the old frame do not change qualitatively 
the attractors found so far. Also, it is not difficult to check that fip* — K^p*/3H^ and are invariant 
under conformal transformation, so that H.^ — fi^ . 

It can be shown that the choice V ~ is the only one that allows scaling attractors in both the old 
and the new metric. Other choices are possible that allow scaling solutions either in the old or in the new 
metric: for instance, V ^ (log/)^^ gives scaling attractors in the new metric but not in the old one. 



V. COSMOLOGICAL PROPERTIES 



Once we have the analytical expression of the attractors, we must consider whether they are viable as 
cosmological solutions. The attractor solution a is inflationary (accelerated) and the scalar field is asymp- 
totically the dominating component. As such, it may match the observations of an accelerated expansion; 
for instance, the value « 0.4 suggested in Ref. Q implies 

M K 0.7 or 3.3. 



On the other hand, since U,^ ^ 1, in order to allow for a substantial fraction in the ordinary matter 
component at the present, the attractor has not to be already reached. 

The solution b has some drawbacks. First, is not accelerated at all for w — 1; second, the constraints 
from nucleosynthesis do not allow a large fraction of energy density in the scalar field, so that it cannot 
provide closure energy. However, as argued in [||, models which reach this attractor compare favorably with 
observations of large scale structure, and may have a simple interpretation in terms of fundamental physics. 

Both solutions a and b are heavily constrained by the upper limits on the variability of the gravitational 
constant. We have 



|G/G| = I///I = 



2 1 



(54) 



Comparing with the observational constraint IG/Gj < alO ^^yr ^, and assuming i w 10 Gyr, we obtain the 
condition 



H > 



1. 



(55) 



Current constraints (see e.g. ) give a « 0.1 or smaller. This implies /i > 20, too large for the attractor 
a to exist. A similar problem arises if / is quadratic. Along the attractor 6, the energy density in the scalar 
field is a constant fraction of the total energy. In MDE (and for /3 = 4 — 3w — 1) this is 



19 + 2^ 



(1 + 2/i) 



2 • 



(56) 



The constraint ( |35|) gives 



< 0.035, 



(57) 



which confines the scalar field contribution to that of a minor component. This constraint is three or four 
times stronger than that imposed by the nucleosynthesis |3] on a minimally coupled field. 



VI. CONCLUSIONS 



In this paper we have investigated a large class of NMC models in the limit of strong coupling with a 
perfect fluid matter component, searching for attractors that might provide a decaying cosmological constant. 
These models include all the cases in which the potential V{4>) is a power of the coupling function /(</>), 
regardless of their functional form. We have shown that 



1. The NMC system can be reduced to a scalar field with an exponential potential, a minimal coupling 
to gravity, and an extra coupling to the matter. 



2. For each pair of the parameters there is one out of four possible scahng attractors: one, a, scalar 
field dominated and possibly accelerated; one, b, decelerated ii w > 0.91 and with constant ratio of 
scalar field total energy to matter; one, c, always decelerated and with constant ratio of scalar field 
kinetic energy to matter; and finally one, d, also always decelerated, and dominated by the field kinetic 
energy. 

3. Attractors c and d are acceptable only in the rescaled frame; in the original frame they lead to a 
gravitational constant either vanishing or infinite. 

4. This choice V ^ is the only choice (in the strong coupling regime) for which there is a scaling 
attractor both in the original and in the rescaled metric. 

5. The constraint on the time variability of G rules out the accelerated models, and only allows a very 
small fraction of the energy density to be in the NMC scalar field. 

Clearly, this analysis is not yet conclusive. Viable solutions might exist for which one or more of the 
following is true: a) the attractors are not yet reached; b) V does not equal /*^; c) the strong coupling 
regime does not apply. For instance, assuming / = 1 + k^^c/)"^, and in the limit of weak coupling , the \G/G\ 
bound can be satisfied for small ^, and the solutions are cosmologically acceptable, although by construction 
do not add much to the minimally coupled model. 
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VIII. FIGURE CAPTION 



Fig. 1. 

Regions of existence and stability in the parameter space w, ^ . In all panels, the eigenvalues of the 
linearization matrix change sign across the thick lines. The color code is as follows: white regions indicate 
that the critical point does not exists; light gray regions, the point is unstable; dark gray regions, the point 
is stable. The dotted lines are at w = 4/3 and w — 5/3 and are useful landmarks in the parameter space. 
From top to bottom, parameter spaces of the critical points a, b, c and d. 

Fig. 2. 

Regions of stability in the parameter spaces. Each region is labelled by the critical point that is stable in 
that region. The gray area indicates where the attractor is accelerated. The stars mark the values of the 
parameters for which we display in the next plot the phase space. 

Fig. 3. 

Phase spaces for various values of the parameters, corresponding to the points marked as stars in Fig. 2. 
The phase space is contained in the positive unitary semicircle. While the phase space of attractors a and 
b are qualitatively similar to those displayed in Copeland et al. |10|, the phase space of attractors c and d 
have no counterpart for (3 = 0. 
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